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Photonic quantum simulators are promising candidates for providing insight into other small- to 
medium-sized quantum systems. Recent experiments have shown that photonic quantum systems have the 
advantage to exploit quantum interference for the quantum simulation of the ground state of Heisenberg 
spin systems. Here we experimentally characterize this quantum interference at a tuneable beam splitter and 
further investigate the measurement-induced interactions of a simulated four-spin system by comparing the 
entanglement dynamics using pairwise concurrence. We also study theoretically a four- site square lattice 
with next-nearest neighbor interactions and a six-site checkerboard lattice, which might be in reach of 
current technology. 

More than a quarter of a century ago, Richard Feynman 1,2 envisioned that a well -controlled quantum 
mechanical system can be used for the efficient simulation of other quantum systems and thus is capable 
of calculating properties that are unfeasible for classical computers. Quantum simulation promises 
potential returns in understanding detailed quantum phenomenon of inaccessible quantum systems 3 , from 
molecular structure to the behavior of high-temperature superconductors. Moreover, quantum simulations 
are conjectured to be less demanding than quantum computations by being less stringent on explicit gate 
operations or error correction 4,5 . Due to these reasons, quantum simulation has led to many theoretical propo- 
sals 6-9 . Recently, various quantum simulators based on different physical platform are being constructed, such as 
atoms in optical lattice 10 " 15 , trapped ions 16 " 21 , nuclear magnetic resonance 22 ' 23 , superconducting circuits 24 , as well 
as single photons 25 " 33 . For recent reviews, see ref. 34-37. 

Motivated by the seminal work of Knill, Laflamme and Milburn 38 photons are proven to be a suitable system for 
efficient quantum computing and quantum simulation. Precise single-photon manipulations and tunable mea- 
surement-induced two-photon interactions are the essential ingredients for photonic analog quantum simulation 
and have been demonstrated. In addition to the high-level quantum control, such photonic quantum simulators 
can produce exotic entangled states which are important for understanding the many-body dynamics in quantum 
chemistry and solid-state physics 8 ' 9 ' 31 . Note that the recently proposed universal quantum computation based on 
coherent photon conversion provides an alternative avenue for achieving photon-photon interaction 39 . 

In ref. 31, a spin tetramer, the simplest example of the Heisenberg-interacting spin system, has been investi- 
gated. In this work we will not only present the detailed study on the tunable measurement-induced photon- 
photon interaction, but also further investigate more complicated Heisenberg models for the four-site square 
lattice with nearest and next-nearest couplings, and the six-site checkerboard lattice. 

Analog quantum simulations with photons and linear optics 

In the case of the simulation of spin- 1/2 particles the photon's polarization is ideally suited as horizontally- 
polarized states \H) and vertically-polarized states \V) represent for example spin-up and spin-down states. 

Moreover, the ability to prepare symmetric polarization -entangled states \^>)^ 2 = ~~r (l-^)i 1-^)2 ± I^Oil^)' 

V2 
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\ x ¥)i2= -j={\H) 1 \V) 2 + \V) 1 \H) 2 ) and the anti- symmetric state 
v2 

\ x ¥-) l2 = ^=(\H) l \V) 2 -\V) l \H) 2 ) enables the establishment of 
V2 

states with bosonic and fermionic character 30,31,40 . The latter shares 
the same quantum correlations as Heisenberg-interacting spins or 
so-called valence bond states 41 . 

The theoretical investigation of strongly-correlated spin systems 
has led to few exact theorems which in some cases are of importance 
for the quantum simulation of chemical and physical models. In the 
particular case of a nearest-neighbor antiferromagnetic Heisenberg- 
interacting spin system it was shown by Marshall 42 that the ground 
state for N spins on a bipartite lattice has total spin zero. This the- 
orem and its extension 43 lead to the fact that the ground state must be 
built as a linear superposition of singlet spin states or valence bonds. 
This constraint, that forces the ground state's total spin to be zero, 
gives rise to various valence-bond configurations that are either 
localized or fluctuating as superposition of different singlet partition- 
ings. Localized configurations are typically referred to valence bond 
solids and delocalized valence-bond states correspond to frustrated 
quantum spin liquids or resonating valence-bond states 44,45 . 
Recently, the photonic quantum simulation of a four spin- 1/2 square 
lattice as an archetype system 31 showed that quantum monogamy 46,47 
plays an important role in frustrated Heisenberg spin systems. 

The experimental setup for studying variable measurement- 
induced interactions is shown in Fig. 1. Our pump source is a 
mode-locked Ti:sapphire femto-second laser with a pulse duration 
of 140 fs and a repetition rate of 80 MHz. The central wavelength of 
the pump is at 808 nm. Then we use a ^-barium borate crystal 
(BBOO) to up-convert the pump pulses to ultraviolet (UV) pulses 



via second harmonic generation. The up-converted UV pulses' cent- 
ral wavelength is 404 nm with a pump power of 700 mW. Then we 
clean the UV pulses with several dichroic mirrors (DM). Photons 1 
and 2 are generated from a BBO crystal (BBOl) via spontaneous 
parametric down conversion (SPDC) in a non-collinear type-II 
phase matching configuration and in an polarization-entangled state 
after walk- off compensation 48 . Photons 3 and 4 are generated from 
BB02 in a collinear type-II phase matching configuration and are 
separated by a polarizing beam splitter (PBS). We guide photons 1 
and 3 to a tunable directional coupler (TDC) and then detect them by 
avalanche photodiodes (APD), which together enable the tunable 
measurement-induced photon-photon interactions to happen. The 
relative temporal delay between the photons is adjusted with a 
motorized translating stage mounted on the fiber coupler of photon 
1. Fiber polarization controllers (PC) are employed to eliminate the 
polarization distinguishability of the two interfering photons. 

Bunching due to the two-photon Hong-Ou-Mandel (HOM) inter- 
ference 49 and the corresponding anti-bunching effect 40 are crucial for 
many quantum information processing protocols, especially for 
photonic quantum computation experiments (C-phase gate 50,51 , 
entanglement swapping 52,53 , etc.), as well as for our photonic 
quantum simulation 31 . The bosonic nature of the photons shows 
up when the input states are superimposed at the TDC such that a 
detection, even in principle, cannot distinguish either of them. This 
leads to a superposition of double occupations on both outgoing 
modes and thus suppression of the coincidence detections, where 
one photon is detected in each output mode. The visibility of this 
HOM dip is one when the TDC is set to have equal splitting of 
transmitted and reflected photons similar to a 50/50 beam splitter. 
As soon as the two input photons can be partially distinguished by 
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Figure 1 | Experimental setup for testing variable measurement-induced interaction between two independent photons. Ultraviolet (UV) femtosecond 
laser pulses are up-converted from a mode-locked Tiisapphire laser. The laser pulses are up converted (BBOO) and the obtained temporal duration and 
repetition rate is 140 fs and 80 MHz, respectively. The UV pulses and the remaining fundamental pulses are separated with several dichroic mirrors 
(DM). Lenses are used to focus and collimate the UV pulses in order to achieve good up and down conversion efficiencies. One pair of the polarization 
entangled photons (photons 1 and 2) is generated from a ^-barium borate crystal (BBO 1 ) via spontaneous parametric down conversion (SPDC) in a non- 
collinear type-II phase matching configuration. Two half-wave plates (HWP) and compensating BBO crystals are used to counter walk-off effects in the 
down-conversion crystal. Another pair of correlated photons (photons 3 and 4) is generated from BB02 in a collinear type-II phase matching 
configuration and are separated by a polarizing beam splitter (PBS). Single-mode fibers and interference filters (IF) are used to clean their spatial and 
spectral modes, respectively. We vary the path length difference between two photons with a motorized translation stage mounted on the fiber coupling 
stage of photon 1. Fiber polarization controllers (PC) are employed to eliminate the polarization distinguishability of the two interfering photons. The 
tunable measurement-induced interaction among photons is achieved by a tunable directional coupler (TDC), followed by a projective measurement of 
one photon in each of the four output modes. 
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unbalancing the splitting ratio the visibility decreases. The depend- 
ence of the ideal visibility (V idea i) upon the reflectivity of TDC (77) is 
the following: 

2*7(1-77) 

Videai = Z Z ~~Z ? • (1) 

l—2rj-\-2r] 2 

In Fig. 2 this reflectivity dependent ideal visibility is plotted in the 
black solid curve. Experimental imperfections due to high- order 
emissions from SPDC and group velocity mismatch reduce Videai- 
The measured visibilities (black squares) are also shown as well as the 
corresponding fit (dashed red curve). 

Pair-wise entanglement dynamics. The main advantage of the 
precise quantum control of individual particles is that inter-particle 
entanglement dynamics can be investigated. By using a similar 
experimental configuration as in Ref. 31, we study the entangle- 
ment distributions among different particles with respect to the 



effective interaction strength that was tuned by the TDC. For the 
quantification of the bipartite entanglement in our system, we use 
the measure of concurrence 54 , which, for a given state p, is 

C(p) = max^0,\/ r Ji — \fhi — \fh> — V^!' where X { are the eigen- 
values of the matrix pHp T H in non- increasing order by magnitude 
with S = cry(x)cry, where a Y = — z|0)(l| + /|1)(0|. While the previous 
quantum simulation characterized the pair-wise energy dynamics of 
Heisenberg interactions that were directly extractable from measured 
coincidence counts, this experiment requires the reconstruction of the 
density matrices to obtain the concurrence values. In the experiment 
we tune the reflectivity of the TDC and hence vary the photon-photon 
interaction strength. Various four-photon quantum states are 
tomographically measured and the density matrices of them are 
reconstructed 55,56 . The concurrence of the two-photon subsystems is 
calculated from the four-photon density matrices by tracing out the 
other two photons. 
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Figure 2 | Experimental demonstration of variable measurement- induced interaction between two independent photons, (a), (b) and (c) Measured 
four-fold coincidence counts are plotted versus the relative optical delay between the interfering photons and fitted with a Gaussian function, when the 
reflectivities of the TDC are tuned to be 0.17, 0.5 and 0.67, respectively, (d) Visibility of the Hong-Ou-Mandel (HOM) dip when measuring four-fold 
coincidences using a TDC. The black squares are the experimental results. The black solid curve is a theoretical prediction based on Eq. 1 and the red 
dashed curve is a fit with the only free parameter of V 5y5 , which is about 0.853. The main reason of this non-ideal value is due to the higher-order emission 
generated from SPDC sources. The error bars are based on a Poissonian statistics. 
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Figure 3 | Measured concurrence for various ground state 
configurations. Entanglement distribution between photons 1 & 2 (black 
square, fitted with black curve), 1 & 3 (red circle, fitted with red curve) and 
1 & 4 (blue triangle, fitted with blue curve) are shown as a function of the 
TDC parameter 6. The rise of entanglement in one pair causes the fall of the 
entanglement in other pairs according to monogamy. The monogamy 
gives rise of the entanglement sudden death and sudden birth 57 " 59 . We 
derive the uncertainties in concurrence from the density matrices, which 
are calculated using a Monte Carlo routine and assumed Poissonian errors. 

Due to the quantum monogamy relations 31 ' 46 ' 47 the total amount of 
pair-wise entanglement stays constant while the change of inter- 
action strength affects the distribution and thus the ground state 
configurations. When tracking the change in entanglement by using 
concurrence, the expected "sudden death" and "sudden birth" of 
entanglement 58 can be seen (Fig. 3). While this concept is typically 
used for studying environment- induced decoherence 59 , a similar 
behavior can be observed here too. In fact our tunable interactions 
allow to mimic a controlled interaction with the environment of two 
additional particles, which opens the possibility to obtain insights 
into complex decoherence mechanisms. We explicitly show the con- 
currence of different photon -pair configuations with respect to the 
TDC angle 0. The relation of the TDC angle and its reflectivity is 
given by 9 = arctan y/rj. One can see that the concurrence for one 
photon pair (e.g. C i4 ) decreases rapidly as we increase 9. At 0 = 
0.274, all of a sudden the entanglement between another spin pair 
(e.g. C 13 ) is born at the cost of the reduced concurrence C 14 , which 
vanishes as 9 is further increased. The observation of similar dis- 
appearance or emergence of entanglement among the other photons 
demonstrates the capability of our quantum simulator for manip- 
ulating the quantum correlations 57 " 59 . 



Figure 4 | (a) A four-site square lattice with / x being horizontal coupling, J 2 
being vertical coupling and / 3 being diagonal coupling, (b) A six-site 
checkerboard lattice. ] x is the coupling between the nearest neighbors and 
J 2 is the coupling on the cross bonds. 

Generalized heisenberg spin model on a four-site square lattice 
and a six-site checkerboard lattice. For the four- site square lattice 
we extend the model in Ref. 31 by adding a next-nearest neighbor 
interaction term (Fig. 4a). The Hamiltonian for this system is: 

H=/ 1 (SrS 2 + S3-S4)+/2(S 1 -S3 + S 2 -S4)+/3(SrS4 + S 2 -S 3 ) 

where S,- is the Pauli spin operator for spin i and J u J 2 and / 3 are the 
coupling strength parameters between different spins, respectively 
(see Fig. 4a). 

We consider the antiferromagnetic case with the couplings / l5 / 2 , 
/ 3 > 0. For each of the three terms in H, the ground state is a pair 
of singlets, |0 = ) = |^-)|¥->, |<D,|) = |^f 3 )|^) and |<D X ) = 
| x F[ 4 )| x F 2 ^) = |O = )-|0)||) (up to normalization). In Ref. 31 we 
showed that by tuning J 2 IJ\ from 0 to °° 5 the ground state, 0^, gradu- 
ally changes from |0=) to |0||). 

With the introduction of / 3 , the ground state of the system is still a 
superposition |0^) = a|0=) + /?|0||) with normalization condition 

2(|a| 2 +|/?| 2 +|a + /?| 2 ) = l. 

Remarkably, tuning J 3 IJi can induce sharp phase transitions with a 
sudden change of the ground state configuration due to the compet- 
ing of the valence bond configurations. In Fig. 5, the ground state 
configurations for different coupling regimes are shown. There are 
three particular interesting phase transitions: 

• For/! = / 3 , | ® g ) suddenly changes from |0=) + |<Xx<> to |0||)when 
J 2 /Ji is tuned across 1 

• For J 2 = J 3 , | %) suddenly changes from |0=)to|0||)-|0><) when 
J 2 /Ji is tuned across 1 

• For/! = / 2 , 1 0^) suddenly changes from |0 = ) + |0||)to |0><)when 
/ 3 //i is tuned across 1 

The case with h = J 2 ^ / 3 is widely studied for square lattice 
systems due to its relevance to cuprates, Fe-based superconductors, 
and other materials 60 " 63 . Previous studies have shown that in the 
thermodynamic limit, when / 3 > / l5 the system is in a diagonal 
Neel ordered state 64 , which is consistent with the ground state |0<) 
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Figure 5 | Ground state phase diagram for various parameter regimes. The coefficients of the ground state, |O g ) = a|0=) + /?|0||)> are shown for 
different ratios of J 2 IJ\ and /3//1. In (a), (b) and (c), we show |a|, \ft\ and |a| + \/3\ y respectively. A clear phase transition can be seen along the dashed line a) 
for h = / 3 , b) for J 2 = / 3 , and c) for ] x = J 2 . 
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Figure 6 | Energy spectrum of the low lying eigenstates of the six-site 
checkerboard lattice system. The energy levels are shown as a function of 
the ratio J 2 IJ\. For the ratio J 2 I]\ = 1 an avoid level crossing of the ground 
state occurs. 

for a minimum of four sites, as discussed above. In the regime where 
J 2 < Ji the configuration for the ground state is still under debate due 
to regions that appear to be non magnetic. Numerical calculation 
have recently shown that this region is highly likely to be a quantum 
spin liquid with Z 2 topological order 64 . However, for only four spins, 
the region with J 2 < J\ has only a single ground state configuration 
|<D=> + |<D,|>. 

We also investigate theoretically the ground states of a ] x — J 2 
Heisenberg model on a six-site checkerboard lattice. The geometry 
of this system is shown in Fig. 4b, where ] x is the coupling strength 
between the nearest neighbor sites and J 2 is the coupling strength on 
the cross bonds. The coupling ratio / 2 //i is the only parameter of this 
system. The introduction of the next nearest neighbor coupling J 2 
makes this system a simple frustrated magnetic model. The study of 
this model on thermodynamics limits is motivated by the three- 
dimensional pyrochlore materials 65 ' 66 . The two-dimensional model 
has been studied by several groups 67 " 73 . It is known that for the regime 
where the coupling ratio hjh^X the system has a co-linear Neel 
order. At / 2 //i ~ 1 numerical calculation 68 suggests a plaquette 



a) 




Figure 7 | The independent ground state configurations of the six- site 

(b)|i^ 2 )> (c)|i^ 3 ), and (d)|i^ 4 ) in our system. 



valence bond solid ground state while the ground state at JilJ\^>\ 
remains under debate 67 ' 69,70 . For this reason future photonic quantum 
simulations might provide answers to these open questions. 

Numerical studies of the six- site checkerboard system were done 
by diagonalizing the Hamiltonian with open boundary conditions. In 
Fig. 6 the energy spectrum of the six low lying states as a function of 
the ratio / 2 //i in the ^^.5^ = 0 subspace is presented. The energy 
spectrum shows an avoid level crossing around / 2 //i ~ 1, which 
indicates a dramatic change of the ground state properties there. In 
analogy to the ] x — J 2 Heisenberg model on a plaquette, where the 
ground states can be expressed as linear superposition of two differ- 
ent dimer coverings whose coefficients depend on the ratio J 2 /Ji, the 
ground states for the six- site checkerboard lattice can also be 
expressed as superpositions of various dimer configurations. For 
the discussed six-site lattice system, fifteen different dimer coverings 
exist. However, only six are independent. By taking the symmetry of 
our system into consideration only four out of the six coverings are 
allowed (Fig. 7). 

Therefore, the ground states can be described as superpositions of 
four dimer coverings, 1^), |i/^ 2 ), |i^ 3 ), |i^ 4 ) in all region of J 2 /J v In 
Fig. 8 we show the contribution of each dimer configuration with 
respect to the coupling ratio J 2 /J v At J 2 /J x = 1 the coefficients for x ) 
and \ \l/ 3 ) are equal and coefficients for \\jj 2 ) and |i^ 4 ) are exact zero. In 
our current convention of the dimer covering wave function this 
particular superposition gives us a plaquette state on the right four 
sites as it is shown in the inset. This coincides with the plaquette 
valence bond solid state in an infinite system. Remarkably, the six- 
site checkerboard lattice thus provides already a valuable hint for the 
true ground state in the thermodynamic limit. The ground states at 
the ratio J 2 > J x have close to equal contribution from and \\j/ 2 ), 
which suggests a significant contribution from a cross-dimer state on 
the left plaquette. We would like to mention that a cross-dimer 
ground state has also been suggested as a potential ground state for 
large J 2 couplings 69 . It is interesting that only a small contribution 
from the dimer configuration \i// 4 ) can be found in all the possible 
cases. 

Summary and outlook. In conclusion, today's available photonic 
quantum technology is reaching the stage where significant advan- 
tages arise for the simulation of particular interesting questions in 
solid-state physics and quantum chemistry. Therefore photonic 
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lattice. The four different dimer coverings are labeled as (a)|^x), 
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Figure 8 | The six-site checkerboard lattice ground state. The according 
superposition of the four independent dimer configurations, li/^), li// 2 )> 
I^X ^4) is shown as function of the ratio J 2 IJ\. The inset shows the 
plaquette ground state at J 2 IJ\ — 1. 

quantum simulations provide exciting opportunities to cover, for 
example, the direct construction of custom-tailored many-body 
wave functions. Impressively, the usage of optical elements such 
as tunable non-polarizing or polarizing beam splitters enables 
entangled few-photon states to construct many-body valence bond 
wave functions for molecular and solid-state systems due to non- 
classical interferences. As we have shown above, this is of particular 
interest in condensed matter physics as it provides insight into the 
frustration of strongly- correlated spin systems and the onset of 
quantum phase transitions. On the other hand, in quantum 
chemistry it effectively allows studying delocalized bonds in 
chemical structures and chemical reactions 6 . Thus being able to 
monitor the full dynamics of individual particles and bonds 
provide some fascinating perspective for the quantum simulation 
of small molecules or reactive centers. 

The main future challenge will be to increase the number of 
photons or degrees of freedom to realize a sufficient amount of qubits 
such that quantum computers can outperform their classical counter 
parts. In general, up to the level of approximately twenty qubits it 
presently appears possible to conceive a system based on bare phys- 
ical qubits. However, given the current experimental limitations, 
operational fidelities and noise sources, it seems that useful system 
consisting of more than twenty qubits could not be realized without 
some level of error correction. But in contrast to the implementations 
of well-known quantum algorithms, such as Shor's algorithm for 
a computationally relevant key-length, the requirements for fault- 
tolerance are much less demanding. As such an open problem relates 
to how much error correction is needed to achieve a useful quantum 
simulation. 

Recent work 6 ' 27 has shown that quantum systems with less than a 
dozen physical qubits are capable of simulating chemical systems 
with a precision that cannot be achieved by conventional computers, 
when processed via almost a thousand discrete gate operations. 
Although such small quantum systems are feasible by using present 
quantum technology, the requirements in terms of gate operations is 
tremendous. Thus, analog quantum gate operations look promising 
in reducing the technical complexity of performing such quantum 
simulation experiments by requiring fewer number of physical gates. 
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